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Abstract. We exploit the gauge-invariant formalism to analyse the perturba- 
tive behaviour of two cosmological models based on the generalized Chaplygin gas 
describing both dark matter and dark energy in the present Universe. In the first 
model we consider the generalized Chaplygin gas alone, while in the second one 
we add a baryon component to it. We extend our analysis also into the parameter 
range a > 1, where the generalized Chaplygin gas sound velocity can be larger 
than that of light. 

In the first model we find that the matter power spectrum is compatible with 
the observed one only for a < 10~^, which makes the generalized Chaplygin gas 
practically indistinguishable from ACDM. 

In the second model we study the evolution of inhomogeneities of the baryon 
component. The theoretical power spectrum is in good agreement with the 
observed one for almost all values of a. However, the growth of inhomogeneities 
seems to be particularly favoured either for sufficiently small values of a. or for 
a > 3. Thus, it appears that the viability of the generalized Chaplygin gas as a 
cosmological model is stronger when its sound velocity is superluminal. We show 
that in this case the generalized Chaplygin gas equation of state can be changed in 
an unobservable region in such a way that its equivalent fc-essence microscopical 
model has no problems with causality. 
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1. Introduction 

Together with quintessence (a scalar field with some potential minimally coupled to 
gravity) [H [21 [3l [4] and other physical and geometrical (modified gravity) models 
of dark energy (see, e.g., the reviews [H [5] for exact definitions), the generalized 
Chaplygin gas [7] (hereafter gCg) is one widely studied model among those proposed 
to describe the observed accelerated expansion of the universe [H [9] . In its generalized 
form, the Chaplygin gas is a barotropic fluid with the following equation of state: 

Pch = -Apc^, (1) 

where A and a are positive constants (for the original Chaplygin gas a = 1). The 
generalized equation of state H]) has been introduced in [7J and analysed in (TUl HH [H] 
and in many other subsequent papers. In contrast to many models describing dark 
energy alone, the gCg gives a unified description of dark matter and dark energy, 
enrolling itself in the class of so-called quartessence or unified dark matter (UDM) 
cosmological models. It allows to interpolate between a dust-dominated phase of the 
evolution of the Universe in the past and an accelerated one at recent time, see Eq. 
^ below. That is why the gCg has attracted much attention in cosmology. 

The evolution of the energy density as a function of the cosmic scale factor a{t) 
of the Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmological model is easily 
obtained from |T]) and from the energy conservation equation: 

^ch=(A+-^y", (2) 

where B is an integration constant usually chosen to be positive (if B < 0, the weak 
energy condition is violated and phantom cosmology takes place, see for example 

mm)- 

From ^ and the parameter wch = Pch/pch and the square sound velocity 
Cg(^[j have the following expressions: 

(3 a) 



■-sCh 



- 1 

dpch 
dpch 

,2 



B 1 



: ~awch- (3&) 



Note that for a > 1, c^qj^ may exceed the velocity of light (set to unity in our notations) 
in the course of the recent or future evolution of the Universe, while it was non- 
relativistic at large redshifts, during the matter-dominated stage, for all values of 
a. 

Cosmological models based either on one or on more than one fluid are essentially 
different. The reason for this is that the various fluid components interact indirectly 
via geometry and may also interact by direct energy exchange. Direct non- 
gravitational interactions are typically neglected but the indirect ones are always 
present and can give rise to totally different evolutions, especially at the perturbative 
level. We will see how this works for the gCg together with baryons. 

The pure gCg, namely the cosmological model based on the gCg alone, has passed 
many tests of standard cosmology. Quite recently, updated constraints for the gCg 
parameters have been pubhshed [HI [TH [17] . However, the behaviour of the gCg under 
perturbations is still problematic [HI [19] . In [19] the study of the power spectrum of 
large scale structures seems to indicate that the best fit value of a is very close to zero. 
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rendering the gCg indistinguishable from ACDM [the latter is indeed the limit of gCg 
for a ^ 0, as one can infer from (|^D|) below]. This feature has to be attributed to the 
gCg sound velocity which, during the cosmic evolution, grows from to ^/a driving 
inhomogeneities to oscillate (if a > 0) or to blow-up (if a < 0). This characteristic 
seems to be common to all UDM models |19) which, as a consequence, might appear 
to be ruled out. 

However, as shown in [20| . it appears that adding a baryon component to the 
gCg the problem disappears and the model is in agreement with observations. But 
the debate still goes on and the authors of [T^ have corroborated their claim about 
ruling out gCg and, in general, all quartessence models through gravitational lensing 
measurements on the basis of the current value of the cosmological parameter as (the 
rms mass dispersion on a sphere of radius 8h~^ Mpc). 

Still, the parameter dg need not necessarily be a confident discriminant among 
different cosmological models since it is clearly a quantity which is strongly affected 
by non-linear effects, while all the calculations performed in [THj are carried out in 
the linear regime. This point was also raised much earlier in [21j in a more general 
context. 

In the present paper we investigate the issue of cosmological perturbations in gCg 
models using the gauge-invariant formalism [22] . We consider both the gCg alone and 
in the presence of baryons, confirming the results of [12] and [50] • In addition, the 
gCg seems to favour structure formation when its sound velocity is superluminal, in 
particular when a>3. 

The structure of the paper is the following: 

(i) In section [5] we briefly outline Bardeen gauge-invariant formalism and derive the 
equations that we will numerically solve in the rest of the paper. 

(ii) In section [3] we investigate structure formation in a gCg-dominated universe and 
display the role of the sound velocity. 

(iii) In section [4] we calculate the theoretical power spectrum in a gCg-dominated 
universe and compare it with the observed one. 

(iv) In section [5] we consider a two-fluid model based on gCg in presence of baryons. 
For this variant of the model we perform the same analysis of the previous two 
sections and compare the results. 

(v) In section [S] we present our conclusions and discuss in more detail why a possible 
superluminal sound velocity of the gCg is not prohibited by causality arguments. 

(vi) For completeness, we add an appendix in which, using the technique developed 
in [23j we write down the exact solution for perturbations of dust-like matter 
and generalized Chaplygin gas in the Newtonian approximation at the matter- 
dominated stage. 

2. The Gauge Invariant Formalism 

The cosmological perturbations issue was first tackled and studied by Lifshitz in 1946 
in the synchronous gauge [24] (see also [25]). In 1980 Bardeen developed gauge 
invariant perturbation theory [22j by constructing suitable combinations of metric 
and stress-energy tensor perturbations which are invariant under a generic gauge 
transformation. The gauge-invariant formalism has been used and reviewed by many 
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authors. Here we follow the notation of 26J and study only scalar perturbations. The 
first-order Einstein equations have the following form: 

A* - Sn (n<P + + S/C^- = a^5p 
n<P + '^' ^ a{p + p)V 

+ n^' + {2n' + $ + 2?t:*' - /c* + 1 a ($ - sy W 

where a{ri) is the scale factor as a function of the conformal time 77. Its present value is 
normalized to unity; 7i(??) = a'/a, where the prime denotes derivation with respect to 
the conformal time; ^{x,r]) and ^'(03,77) are the Bardeen gauge-invariant potentials; 
6p{x,ri) and 6p{x,ri) are the gauge-invariant expressions of the perturbations of the 
energy density and pressure; V(x, rf) is the gauge-invariant expression of the scalar 
potential of the velocity field; a{x,r)) represents the shear. Finally, 7^ {i,j — 1,2,3) 
is the spatial part of the FLRW metric and K. — 0, ±1 is its curvature parameter which 
corresponds respectively to fiat, close and open geometry. The vertical bar denotes 
covariant derivation with respect to jij. Units are chosen such that AttG — c— 1. See 
[26] for a detailed derivation of system (j4|) . 

We introduce two assumptions which simplify (jj]): 

• We assume /C = 0; 

• We neglect shear perturbations, namely we assume cr = 0. This implies $ = 

For a cosmological model based on N non-interacting fluids we write for the 
background quantities: 

while for the perturbations in the linear regime: 

= Ejll ^P»> = Ejll ^P^■ (6) 
The mixed time-space component of the perturbed energy-momentum tensor will 
look like X^iLi (Pi Pi) ^1 since each fluid will contribute with its own velocity. 

Adiabatic perturbations of a barotropic fluid are characterized by the following 
equation of state: 

Spi = cl^Spr, (7) 

where 

2 - ^Pi (Q\ 
Opi 

is the sound velocity at constant entropy for the generic component i. In the 
cosmological perturbation theory it is customary to deal with the density contrast: 

S. ^ ^, (9) 
Pi 

where pi is the background energy density of the component i. Since we work in the 
linear regime it is useful to take the spatial Fourier transform of ([4]) which allows 
us to treat each mode independently. With the above assumptions and trading the 
conformal time for the scale factor a, system (jlj for a generic multi-fluid model is 
rewritten as follows: 

n'^ + an^ = aj:f^^{p,+p,)V, (10) 
{anf $ + (40^2 + a^nn) $ + (2aHH + hA<^^ J^^^i Cs^P^Si, 
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where the dot denotes derivation with respect to the scale factor a (for the sake of 
simpUcity we adopt the same notation for the Fourier transformed quantities as for 
the original ones). 

If iV = 1, the system is determined and can be solved by eliminating S and by 
extracting the following second order equation for $: 

Once solved, we can use the solution of (fTT|) in the first equation of (fTO| to find the 
solution for S. 

Instead, if > 1 we need 2(iV — 1) additional equations to make pU]) determined. 
These equations should take into account the possible interactions and the energy 
exchanges among the different fluids. Since we have assumed that they are not 
mutually interacting directly, the required additional equations to be added can be 
simply chosen as the energy conservation equations, (5T^.^ = 0, given separately for 
each further component. Namely the continuity equation 6Tq.^ = 0: 

6, + - (4 - w,) 6, - 3$ (1 + w,) + ^{l + w,)V, = (12) 
and the Euler equation ST[^^ = 0: 

[(p. +P^)V,]-+^{p,+p,)V.-'^S.- iEl±Pll<p = 0, (13) 
a n n 

where Wi = Pi/pi- Both (fT2l) and are Fourier transformed, that is why there is 

no I subscript in (fT3)l . 

It is possible to write the combination of (|10p . and (I13p as a system of N 

second order equations, one for each 6i. One needs to derive (HH), then to express 

Vi and Vi in terms of the other J^'s through ^TU^ . (fT5)) and again In the case 

N = 2 the resulting system has the following form: 



f ^1 + 

























# + # ) <52 + If 'ji + ^<52 + j^h = 0, 
where the coefficients have the following expressions: 

a a ^ ' 



(14) 



iV? (34 - 1) + fc2 (3c2i + 1) + 6 _ -una] 
-ian{pi+pi) . ■ ^, (15) 



3n^ (34 ~i) + k^ (34 + 1) + 6 _ nUa) 

Bi=-3anp2{l + wi) -. ^, (16) 

Ci = + 4 (nna + (4 - wi) + 3— (cl^)'- 3 (l + u;i) + 

3H2(2 + 3c,\)(pi+pi) Vsi iM ) 



/c4 + 3 (37^2 + fc2) I^^TYa - 7^2^ 
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m"^ -mna] {pi+pi 



37^2 (c^i - wi) 



fc4 + 3 (3^2 + fc2) (7^?ia - H' 



(17) 



L'i = -3(l + wi)c,^2P2 4 
3^2 (2 + 3c2,) p2 (1 + Wi 



fc2 - 3 (( 



■s2 



W2) (3^2 



fc4 + 3 (37i:2 + /c2) (^7i;7i;a _ 

, . fc2 + 3 

+ 3^2 - 6HHa p2 (1 + m) 



fc2+3 


T^T^a - 7^2 _ 


37^2 (c22 _ ^^,2) 


/fc4 + 3(37i2 + A:2) 


'nka - re 





(18) 



The coefficients A2, B2, C2 and D2 have the same form, but with the interchange 
1 <-> 2 in the subscripts. 

As shown in 23J, at the matter-dominated regime, i.e. when pi <^ pi, P2 'ti P2 
and Csi, Cs2 *C 1, the above coefficients are far simpler and system (jl4p can be rewritten 
in the following quasi-Newtonian form: 



' ^1 + 






< 


ii- 










^2 + 


(»- 





+ Jji^clj^Sl = (^l*^! + /'2'^2) 



(19) 



73- {pi5i + P252) ■ 



This approximation is sufficient to study the formation of all gravitationally bound 
objects (galaxies, in particular) which become nonlinear for z > 1. As pointed out in 
[25] . it is possible to recover the system (ITOl) by considering the limit fc 3> H, too, but 
it is not a necessary constraint in the matter-dominated regime when Csi and Cs2 are 
small. It is interesting that the system (IT^ can be exactly solved for the two-fluid 
model considered in Sec. [5] The solution is exhibited in the Appendix. 



3. Perturbations in a Chaplygin gas dominated universe 

In the gCg dominated universe the Hubble parameter scales as follows: 

where A = A/ {A + B) and Tio = Hq is the Hubble constant. 

In a typical analysis of the gCg model it is customary to take a < 1 since the 
square sound velocity tends to a for a — > c», as can be seen from For a — 1 the 
sound velocity will tend to unity in the far future, at t — s- 00. The new regime that 
we investigate in the present paper is a > 1, thereby making gCg superluminal after 
a sufflciently long but still finite time. 

From p&j) it is straightforward to calculate at which redshift the transition to the 
superluminal gCg would occur. For generic values of A and a: 



A{a-l) 



\-A 



1. (21) 



Given a > 1 and A, it is plausible to expect that some, hopefully observable, 
cosmological effect would take place at that redshift. 

It is crucial to point out that a and A are not independent, but linked by (j20p . 
once we have some constraints about the background expansion. In the present paper 
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we consider Zt^, i.e. the redshift at which the transition to the accelerated phase of 
the expansion takes place. We derive its expression from ([^D|) : 



Ztr 



2A 



3(q + 1) 

- 1- (22) 



.1-^ 

It follows from the most recent SN la observations combined with CMB and 
baryon acoustic oscillations (BAO) data that ztr ~ 0.7 with a rather large uncertainty 
(at least ±0.1 or more at the 2a confidence level), see e.g. Fig. 7 in [27] and also 
[ni EH [29]. For comparison, if only SN la are used, then the Gold+HST dataset 
[50] leads to ~ 0.4, see Fig. 2 in [57] and also [21]) while other datasets (SNLS 
and ESSENCE) alone produce values of zu- ^ 0.7, too (see [17] [27]). However, these 
numbers strongly depend on the present value of the non-relativistic matter density 
r^mo which, in turn, has to be expressed through A in the gCg case (the values given 
above correspond to fimo ~ 0.28). Direct calculation of ztr for UDM models [TB] [28] 
shifts this quantity to even larger redshifts: ztr = 0.8 — 0.9. Actually, in the latter case, 
the authors consider only silence quartessence (namely the effective sound velocity is 
assumed to be zero) and the parameter a is constrained to be a = —0.06 ±0.1. 

Then from (|22p we can write down an explicit relation between A and a: 



A=^i±i!l)-_. (23) 
2+(l + zt,)^(^+") 

Note that the uncertainty a^^^ of the transition redshift propagates on A through the 
following formula: 

^ 6(a + l) 

~ [2(1 + Ztr)-3("+l) + 1]2(1 + Ztr)3"+4''^*-' ^^^^ 

therefore the larger is a, the smaller is the bias on A. 

A widely used cosmological parameter stemming from CMB observation is i?, the 
comoving distance to the last scattering surface scaled to fimOj namely: 



dz 



R = HoVn^ — -, (25) 
Jo J^v-) 

where z\s ~ 1089 is the last scattering surface redshift. In UDM models, the quantity 
r^nio should be taken from the asymptotic behaviour of Ti^ for z 1 at the matter- 
dominated stage, so flmo = (l - ^)^^^"^^^ (or f^mo = (l - A)^^^"^^^ (1 - fi{,) + if 
baryons are taken into account, too). At the la confidence level, R = 1.71 ± 0.03 (see 
for example [32]). 

An interesting limiting case is the 'super-duper luminal' one: a ^ oo. Then 
H = H/a = const = Hq for z < Ztr and H^{z) = ((1 + z)/(l + ztr)f for larger z. 
In this case, one can obtain an analytic expression for R: 

z^Y 2 2 

Then the above mentioned observational window for R is reached for ztr = 1-0 ± 0.1. 
Combining ([22]) with ([21]) we find that 



Q, _ 1 \ 3(a + l) 

(l + ^tr)-l. (27) 

If a = 3, then Zg = ztr- For a > 3, > Ztr but Zg approaches Zn once more for 
a oo. For a fixed Ztr, the maximal value of Zg is reached for a « 8.182 when 
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(1 + Zs)/{1 + Ztr) ~ 1.048. Thus, for a > 3, the transitfon from subluminal to 
superluminal gCg occurs approximately at the same time as the transition to the 
accelerated phase of expansion of the Universe. 

We now compare the evolution of perturbations in the gCg-dominated universe for 
different values of a. We choose an integration range which starts at the decoupling 
era, namely z ~ 1100 (a « 10~^), and ends up when the first structures, namely 
protogalaxies, were formed (a ~ 0.1, or z « 10), so that the linear approximation 
holds true. 

We numerically solve pT|) choosing [—1,0] as normalized initial conditions for 
$, $ . In our calculations we use ((23|) to properly choose A as a function of a. We 
show the results using ztr = 0.8 ± 0.16 and ztr = 1.0 ± 0.1. 



u 



2.5 



xlO" 



: 0.8 ±0.16 



2 - 



1.5 



0.5 



-0.5 




0.01 



0.02 



0.03 



0.04 0.05 0.06 
Scale factor a 



0.07 



0.08 



0.09 



0.1 



Figure 1. Evolution profiles of Sch for difi'erent values of a and for k = 100 h 
Mpc~^. The transition redshift is ztr = 0.8 it 0.16 
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Figure 2. Same as figure [l] with , 



1.0 ±0.1 



The plots in figures [T] and [2] display the evolution of the density contrast in the 
gCg, which we call i5ch, for different choices of the parameter a and for k = 100 
h Mpc^^, which corresponds to a scale of order 50 kpc, which is typical of a 
protogalaxy. The case a = corresponds to the ACDM model. 

Indeed, here we are first investigating if formation of gravitationally bound objects 
including galaxies takes place at all. That is why we have chosen so small a scale and 
we have stopped at z = 10. The gCg UDM model, like all CDM-like ones, belongs 
to the class of bottom-up, or hierarchical clustering, scenarios of structure formation 
where small scale structures are formed first. The question of large-scale structure in 
this model at scales which remain linear for z < 1, of course, requires integration up 
to z = 0. It is considered in the next section. 

Part of the results found in [THj are already confirmed by this analysis, without 
a deeper study involving the power spectrum. The development of oscillations takes 
place too early for a > 10~^, thus preventing structure formation. 

However, two further interesting features can be extracted from figures [T] and 
[2] First, it seems that there exists a critical value, namely a « 0.1, for which the 
deviation of the gCg from the a = linear growth is maximal. Second, the larger the 
value of a (in the range a > 0.1) the smaller the deviation from the a = behaviour. 
Indeed, when a > 1, i.e. in the supcrluminal regime, the oscillations are absent and 
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the growth seems even to be enhanced. 

Both the above mentioned features can be explained by the behaviour of the gCg 
sound velocity as a function of a, which is displayed in figure [3] for Zf^ = 0.8 ± 0.16 
(the plots for ztr = 1.0 ± 0.1 are very similar). 



z =0.8 ±0.16 

tr 




a = 3 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 

Scale factor a 



Figure 3. The gCg square sound velocity given in l)3fe| l plotted as a function of 
a for different values of a. A is given by II23I I with ztr = 0.8 ± 0.16. 

For a ~ 0.1 the sound velocity becomes non negligible much earlier than at other 
values of a. The range lO"*" ^ a < 3 appears thus to be ruled out since structure 
formation is prevented. However, to understand what happens out of this range it is 
necessary to study the clustering properties of matter at the present time, namely the 
power spectrum. This will be done in the following section. 
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4. Large Scale Power spectrum of Chaplygin gas perturbations 

In the previous section we have shown that in the gCg model galaxy formation is 
possible only when a is very small or a > 3. In this section we try and constrain these 
bounds further by studying the clustering properties of gCg perturbations, namely the 
power spectrum P{k), and comparing it with observation. 

We exploit the SDSS (Sloan Digital Sky Survey) data as analysed in [33J. The 
sample considered by the authors consists of 205,443 galaxies observed before July 
2002. The data consist of RA (Right Ascension), Dec (Declination) and redshift 
z for each galaxy. The redshift is converted to comoving distance by using a flat 
cosmological model with a cosmological constant Ha — 0.7 (it can be shown that our 
results are robust to this assumption). 

The power spectrum is computed at present time, namely at z — 0. The data 
processed in [33] are displayed in table [TJ 



Table 1. Effective k, window functions and measured values of P{k) with relative 
bias factors and standard deviations, [k] = hMpc~^ and [P{k)] = {h^^ Mpc)^ . 
From the SDSS data reported in |33l . We have neglected two rows of data in 
which a is greater than the respective P{k). 



k 


Low k 


High k 


P{k) 


Bias 


fj 


0.01819 


0.01502 


0.02440 


33254.63896 


1 


16745 


24572.79357 


0.02405 


0.01984 


0.03136 


38360.58264 


1 


16674 


13320.48227 


0.02776 


0.02300 


0.03571 


24143.07699 


1 


16613 


10047.21709 


0.03201 


0.02663 


0.04048 


19709.29306 


1 


16532 


7413.82995 


0.03691 


0.03094 


0.04601 


12595.77528 


1 


16428 


5485.84481 


0.04257 


0.03591 


0.05215 


13558.60084 


1 


16294 


4077.68223 


0.04912 


0.04162 


0.05983 


18311.08054 


1 


16124 


2974.15048 


0.05670 


0.04824 


0.06873 


12080.73574 


1 


15910 


2140.39514 


0.06527 


0.05561 


0.07867 


9217.46084 


1 


15647 


1580.01728 


0.07529 


0.06455 


0.08999 


9750.49986 


1 


15317 


1127.95803 


0.08698 


0.07507 


0.10351 


9529.63935 


1 


14906 


818.49019 


0.10037 


0.08670 


0.11898 


6384.82545 


1 


14409 


601.70752 


0.11581 


0.09985 


0.13679 


5294.92081 


1 


13813 


446.65015 


0.13360 


0.11478 


0.15748 


4629.83669 


1 


13109 


335.07950 


0.15412 


0.13181 


0.18139 


3573.92724 


1 


12290 


254.13547 


0.17774 


0.15114 


0.20891 


3393.82814 


1 


11358 


195.23431 


0.20489 


0.17288 


0.24062 


2298.10000 


1 


10320 


152.73819 


0.23592 


0.19690 


0.27653 


1597.14976 


1 


09193 


123.50981 


0.27062 


0.22219 


0.31395 


1105.42903 


1 


08019 


107.20069 


0.30618 


0.23123 


0.34782 


1012.72489 


1 


06912 


110.33632 



We compute the gCg power spectrum for different values of a, by first solving (jll[) . 
thus finding the transfer function and then applying it to the following prior: 

2 



14 



Nk 



ln(l + 2.349) 



2Mq 

where A'^ is a normalization constant and 
k 



1 + 3.89q + {le.lqf + {5A6qf + (e.llqf 



-1/2 



q 



(28) 



(29) 



where fixo is the cold dark matter (CDM) energy density fraction evaluated today, 
and h is the Hubble constant in units of 100 km s^^ Mpc~^. In our calculations we 
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shall set h — 0.7. Since the gCg has the same asymptotic behaviour as the ACDM 
model in the matter-dominated stage, then the cold dark matter fraction is given by 

nxo={l~A)^. (30) 

The prior ([28]) is obtained by applying the BBKS transfer function [34] to the scale 
invariant Harrison-Zel'dovich spectrum. The normalization constant N is computed 
through a best fit of the data for each a. In our calculation, A is constrained by 
using ztr = 0.8 ± 0.16 and ztr = 1.0 ± 0.1. 



z =0.8 ±0.16 

tr 




10-' 

k [h Mpc"'] 



Figure 4. Theoretical power spectra for different values of a compared with the 
measured one. A is given by I I23I I with ztr = 0.8 it 0.16. Note that the plots for 
a = and a = 10"'' are superposed. We have extracted the envelopes of the 
oscillations of the density contrasts for a = 3 and = 5 and have used them 
for the calculations. For the other cases this procedure was unnecessary since no 
oscillations were present. 
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z = 1.0±0.1 

tr 




10-' 

k [h Mpc"'] 

Figure 5. Same as figure|4] with ztr = 1-0 it 0.1. 

In figures [?] and [5] the power spectra have been computed by considering the 
envelope of the oscihations of the density contrast Sch and by using it in the subsequent 
calculations. This procedure serves to clearly compare the diS'erent slopes of the 
calculated power spectra with respect to the observed one. Moreover, it has been 
carried out only for a — 3 and a = 5. For the other ones it was unnecessary as no 
oscillations were present. 

An interesting feature that can be drawn from figures 2] and [5] is that the larger is 
a the more the power spectrum tends to a limiting behaviour which is systematically 
below that of the ACDM one. This happens because the suppression of the gCg 
transfer function relative to that of the ACDM model behaves approximately as fc~^ 
for large values of a after averaging over oscillations (actually, the exponent is slightly 
larger than —1 because of the accelerated expansion at recent times). 

Figures [6] and [7] show what happens if the oscillations in the density contrasts are 
taken into account. The power spectra for a = 3 and a — 5 are highly irregular and 
their slopes are hardly recognizable. 
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z =0.8 ±0.16 

tr 




10-' 

k [h Mpc"'] 



Figure 6. Theoretical power spectra for different values of a compared with the 
measured one. A is given by ||23|I with ztr = 0.8 ± 0.16. Note that the plots for 
a = and a = 10~^ are superposed. 

In conclusion, the pure gCg does not seem to work from a perturbative viewpoint, 
even when it is superluminal. However, a reasonable possibility is to allow the gCg 
to have non-adiabatic perturbations, which is a natural assumption since it is not a 
pressureless fluid. An attempt in this direction has already been performed in |35j 
and [36] (here also a baryon component is considered). However, the results are based 
on the ad hoc assumption of silent perturbations (namely Spch = 0), which seems 
questionable. 
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Figure 7. Same as figure|6] with ztr = 1-0 it 0.1. 



5. Two-fluid model: generalized Chaplygin gas plus baryons 

In this section we consider a two fluid model including baryons together with gCg. 
The Friedmann equation governing this model has the following form: 

where subscript b refers to the baryonic component and subscript indicates that the 
corresponding quantity is evaluated at the present epoch {a ~ 1). 

System (|14p is solved numerically in the same integration range of the previous 
sections and setting ilbo — 0.04. Before proceeding it is necessary to find a relation 
between A and a which ensures the preservation of the correct background expansion, 
e.g. the acceleration redshift ztr ~ 0.6. 

From (|3ip it is possible to work out the following expression: 

^[(l + ztr)-^<°+^^ + l/2]-l/2 ^ »bo .32^ 

{i[(l + 2tr)-3(o+l)_l] +1}"/"+! 2{l-n^oy 



1-A 



(31) 
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which cannot be inverted to give ^ as a function of a. However, since f^bOi is small 
we can use again to that purpose. 

As in the previous section, we choose (j28|) as a prior for the power spectrum, but 
since this time we have the presence of baryons, ilxo has to be chosen differently. We 
exploit then Sugiyama's shape correction [37] : 



xo 



bO 



{1-n 



60 J 



1 



exp 



2hn 



bO 



60 



+ (1 - ^bo) (1 - A) 1^ 



(33) 



which properly takes into account the presence of baryons. 

In the appendix we show an exact solution, in the matter dominated regime, for 
the evolution of perturbations in the model of the present section. 



z =0.8 ±0.16 

tr 




k [h Mpc"'] 



Figure 8. Power spectra of the baryonic component computed for different values 
of a with ztr = 0.8 ± 0.16. 
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Figure 9. Same as figure|8]with ztr = 1.0 ± 0.1. 

In figures [8] and [9] we show the baryonic power spectra computed for different 
values of a compared to the observed one. We have chosen as normalized initial 



conditions the following values: 



— [1,1,1,1]. We have chosen the 



<5b, (5b, Sch, 6ch 

same initial conditions for all the quantities essentially because after the equivalence 
era the gCg behaved as dust and therefore its density contrast grew linearly (as a 
function of the scale factor) . 

Unlike what happened in the previous section, this time there were no oscillations 
in Sb, since the baryonic density contrast is unaffected by the instabilities of the gCg 
one. In fact, the results are different from those inferred from figures [4] and O The 
agreement with the observed power spectrum is better. It is particularly good again 
for very small a but also in the range a > 1 . Therefore the gCg seems to work better 
either when it is practically indistinguishable from ACDM or when its sound velocity 
is superluminal. 

The good agreement of the theoretical power spectrum with the observed one 
can be explained realizing that gCg perturbations are indistinguishable from those 
of baryons at the beginning of the matter-dominated stage. Therefore they really 
act as dark matter and compel baryons perturbations to grow and form structures. 
On the other hand, when the gCg sound velocity increases, the growth of the baryons 
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inhomogeneities is damped, but does not oscillate. This fact can be observed in figures 
[inianddlj Moreover, structures which have already formed are not influenced by the 
gCg sound velocity, so neither is the large scale power spectrum. This was already 
pointed out in [5D]. 



z =0.8 ±0.16 



10^ 




S 10' 



10" 



10 



10"- 
Scale factor a 



10" 



Figure 10. Evolution profiles of 5]^ for different values of a and for k = 100 h 
Mpc~^. Note that the plot for a = and a = 3 are superposed, ztr = 0.8 ±0.16. 

In figures [TO] and [TT] we display the evolution of 5h for different values of a choosing 

again ^b, <5b, (5ch, ^ch = [Ijlil,!] as normalized initial conditions. As in section [3] 

we have chosen fc = 100 /i Mpc^^. Note that the plots corresponding to a = and 
a = 3 are superposed. 

We infer from figures [TOl and [Til that the gCg sound velocity has a damping effect 
on the growth of baryon inhomogeneities. These are not compelled to oscillate, but 
structure formation is delayed, if not prevented. In order to render structure formation 
as similar as possible to that of ACDM, we need either very small values of a or a > 3. 
The latter range for a just leads to Zg > ztr as was shown in section [31 
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Figure 11. Same as figure [TOl with ztr = 1-0 it 0.1. 

6. Conclusions 

From our analysis it emerges that: 

• The generalized Chaplygin gas cosmological model, with no additional fluid 
components, is compatible with structure formation and large scale structure 
clustering properties (encoded in the power spectrum) only for a sufficiently small 
(a < 10^^), in which case it is practically indistinguishable from the ACDM 
model. 

• Adding to the generalized Chaplygin gas model a baryon component we find 
that the growth of the density contrast and the large scale structure clustering 
properties are compatible with observations for all values of a. However very 
small values of a and the case a ^ 3 are favoured. 

Note, however, that the integrated Sachs- Wolfe (ISW) effect (not considered in 
this paper) leads to rather large enhancement of low CMB multipoles C/ in UDM 
models, in contrast to a modest increase of C; in the ACDM model [H], that results 
in the exclusion of the gCg with 0.01 < a < 1 [TH1I351[1S| (see also the recent paper [17] 
for an analytical treatment of the ISW effect for UDM models). It may be interesting 
to investigate if it remains so large for all a > 1. 
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Now, since most of the new results of this paper refer to the superluminal case 
a > 1, it is natural to consider the question if this range is physically admissible in 
more detail. The fact that the group velocity CgCh in media may exceed the light 
velocity in vacuum (unity in our notations) is not, in itself, unphysical. Not only it 
is not prohibited from the theoretical point of view, but this phenomenon has been 
actually observed in laboratories in case of light propagation through dispersive media. 
As was proved long ago (see e.g. [38]), what is really necessary in order not to violate 
causality is that the signal (or, the wavefront) velocity should not exceed 1. 

Then, what is the signal velocity for the gCg? It does not appear possible 
to answer this question at the hydrodynamic level since the gCg phenomenological 
equation of state in the normal (non-phantom) case requires pch 

> so that 

the limit pch — > corresponding to a wavefront expanding into vacuum cannot be 
studied. Therefore, to determine the value of the signal velocity for the gCg, one has 
to use some underlying microscopic field-theoretical model from which the equation 
of state (dl arises in macroscopic hydrodynamics. 

We will use the so-called tachyon representation of the gCg for this purpose where 
the gCg is described by a scalar field (j) (fc-essence) with the Lagrangian density [11] 



C = PCh = -Vo 



X 



1 

: X = 770,^-^'^ < Vo , (34) 



2' 



PCh = "^X^ - C = Vq 



1 - 



. 1+° 

Vr 



A = ^(,1+" . (35) 



Though other microscopic models are also possible (see e.g. discussion in 39J in the 
case a = 1), the present one has the following advantages: a) its solutions are in 
one-to-one correspondence to those of the hydrodynamic gCg if X > 0, and b) it has 
the same equations for perturbations, too. 

Now we are able to consider the limit X corresponding to a wavefront 
expanding into vacuum with no fc-essence in front of it. To prove that the signal 
velocity is equal to unity in this case, it is sufficient to show that the wave equation 
corresponding to (|34p , namely 

+ = 0, (36) 

has plane- wave solutions of the type f{x — t) in flat space-time (since the signal 
velocity is also the infinite momentum limit of the phase velocity of a wave). It is 
straightforward to check that the sufficient condition for this to hold is that both £x 
and Cxx are finite at X = (a similar remark is made in the recent paper |40| . see 
also |41]). In addition, Cx{X = 0) should be non-negative in order for ghosts not to 
be present. 

In the superluminal case a > 1, the Lagrangian density (|34p does not satisfy this 
condition since Cx and Cxx diverge at X —^ 0. However, without going into a more 
detailed analysis of the singularity structure at a null hyper-surface in this case, we 
propose a simple sufficient way to cure the causality problem for the gCg. Namely, 
let us assume that the expression ([34)) changes to 



C = -Vo+X[^^] (37) 
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for X -^Vi where Vi <C Vo and that it smoothly matches (|34p at X ~ Vi . Then it is 
clear that the signal velocity is exactly unity. The Lagrangian density (|57)) describes 
a massless scalar field minimally coupled to gravity and a cosmological constant. 

The correspondingly corrected macroscopic gCg remains the same in the 
observable region X > Vq, or pch + Pch ^ Vq (note that Vq is of the order of the 
present day critical energy density in the Universe) but will become a mixture of a 
cosmological constant and the extremely stiff ideal fluid {p = p) in the far future 
when X will drop to Vi . With the same accuracy, one can say that the corrected gCg 
becomes the usual Chaplygin gas when its density is very close to Vq, much closer 
than the one which occurs presently. Thus, by slightly changing the gCg equation 
of state in the unobservable region, the causality problem is cured. As a result, the 
present gCg superluminal sound velocity is a transient phenomenon which disappears 
at t —* oo. 

Finally, note that in the limit a ^ oo, the gGg realizes the Parker-Raval scenario 
[12] of a sudden transition from the matter-dominated stage with a{t) oc i^/^ to the de 
Sitter expansion with a constant curvature (modulo a small term cx flf,). The original 
idea of Parker and Raval to produce such a scenario using non-perturbative vacuum 
polarization of a light minimally coupled scalar field is questionable, both regarding 
derivation of the effective action and stability of the resulting specific kind of the f{R) 
theory of gravity (here R is the Ricci scalar). But the same scenario based on the 
limiting case of the gCg does not have these problems. This shows once more that the 
superluminal gCg deserves further, more detailed study. 
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Appendix A. Exact solution for perturbations of the gCg+baryons model 
at the matter dominated regime 

At the matter-dominated regime, the system ([T9l) describing gCg and baryons takes 
the following form: 

^ ^ = ^ {pbSh + PchSch) ^ , ^ 

2 (A.l) 

Ch + -^TfpCgSch — {Pb&h + PCh^Ch) , 




where Cg is the gCg square sound velocity given by (Pa|) and (|36|. 

This system can be exactly solved using the technique developed in [^3] for 
cosmological perturbations of two-fluid models in the Newtonian approximation. Let 
us introduce the following variables: 

x^kj-'^a-i'', 7 = -2a-|, (A.2) 
where a is the gCg parameter. It is possible to extract from system (|A.ip a fourth 
order equation for Sch- 

^ ^ 37) ^ 37) 7 
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Sch = 0, (A.3) 



27-1/3 ^ , Iff 1\ 2 



.3 
where A = x • d/dx. 

As found in [23] , the general solution of (jA.Sp can be represented in terms of the 
Meijer G-functions (the generalized hypergeometric functions): 

+ CsGf, {xe'TblSkJ + C,G\l (xe-'^^l-- , (A.4) 

where: 

ai,2 = ^ (1 T VI + 24r!bo) 

bi = -j^, &2 = (A.5) 
^3 = ^, ^'4 = ^, 

the Ci {i = 1, 2, 3, 4) are integration constants and the Meijer G-functions are defined 
as follows [43] : 

' - 2.i4 uu^^na,-s)uUm^,rii-b, + sf ^^-'^ 

where T{s) is the Gamma function and the integration contour 7l circuits all the poles 
of r(l — Uj + s) anti-clockwise and all the poles of r{bj — s) clockwise (all these poles 
lie along the real axis). 

The solution for Sh has the same functional form as (|A.4p . with: 

1 



p,q 



ai.2 = 1 + — [1 T + 24r!choj , (A.7) 
where flcho — ^ ~ ^bo since we assume a spatially flat universe. 
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